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1.0  INTRODUCTION 

Many  naturally  occurring  flows  as  well  as  those  encountered  in  engineering  applications 
are  time  dependent.  In  such  cases  the  velocity  field  responds  to  an  imposed  pressure  gradient 
or  other  external  excitation  in  an  extremely  complex  fashion.  To  illustrate  the  possible 
effects  of  time-varying  conditions  on  fluid  response,  it  is  useful  to  consider  a  class  of  internal 
flows  that  are  called  (equivalently)  pulsating,  pulsatile,  or  oscillating.  Such  flows  are 
characterized  by  temporally  periodic  variations  of  the  imposed  pressure,  which  may  result 
from  (1)  unstable  combustion  processes,  (2)  acoustic  disturbances,  (3)  vibrations  at  flow 
boundaries,  (4)  the  action  of  reciprocating  pumps,  or  (5)  unstable  pressure  regulators  and 
valves.  It  is  known  from  analysis  (Refs,  1  through  4)  and  has  been  demonstrated 
experimentally  (Refs.  5  through  8)  that  flow  pulsations  may  (1)  promote  transition  to 
turbulence,  (2)  alter  the  turbulent  structure  of  the  flow,  (3)  cause  flow  reversals  during  a 
portion  of  a  cycle,  and  (4)  significantly  augment  or,  in  some  cases,  inhibit  heat  transfer.  The 
velocity  field,  moreover,  varies  in  its  response  to  a  periodic  pressure  disturbance  with  the 
peak  velocity  amplitude  lagging  the  peak  pressure  amplitude  by  as  much  as  90  deg  (i.e.,  one- 
fourth  the  period  of  the  disturbance)  on  the  centerline  of  a  tube.  Perhaps  the  most  striking 
feature  of  bounded,  pulsating  flows  is  the  observation  that  the  maximum  fluid  velocity 
frequently  does  not  occur  on  the  centerline  of  the  duct  or  pipe  so  that  (1)  a  greater  portion  of 
the  mass  flow  is  carried  in  the  annular  region  near  the  wall  and  (2)  the  instantaneous  wall 
shear  stress  may  be  significantly  greater  than  in  steady  flows  at  the  same  Reynolds  number. 

It  is  obvious  that,  even  in  this  relatively  simple  case,  transient  fluid  motions  may 
introduce  effects  that  cannot  be  ignored  in  conducting,  or  analyzing  data  from, 
experimental  studies.  The  complexity  of  the  phenomena  is  such,  moreover,  that  attempting 
to  draw  conclusions  on  the  basis  of  quasi-steady  models  is  not  advisable  and  often  is  a 
serious  error.  In  view  of  these  considerations,  it  is  necessary  to  develop  practical 
computational  models  that  will  allow  ihe  determination  of  the  details  of  a  variety  of 
transient  flows  common  in  AEDC  test  facilities  to  serve  as  a  guide  for  both  the  establishment 
of  test  procedures  and  the  interpretation  of  test  data. 

Undoubtedly,  the  most  general  way  to  obtain  reliable  models  of  transient  fluid 
phenomena  will  be  through  the  development  of  finite  difference  or  finite  element  models 
which  solve  the  governing  partial  differential  equations.  Pending  the  development  of  such 
programs  and  as  an  aid  in  their  evaluation  it  is  useful  to  develop  mathematically  exact 
solutions  to  those  specific  problems  which  lend  themselves  to  exact  analysis.  Various 
analytical  solutions  have  been  proposed  for  the  velocity  field  on  pulsating  flows.  Sexl  (Ref. 
9)  obtained  a  solution  for  the  velocity  field  in  a  tube  in  which  only  a  purely  oscillatory  flow 
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existed.  His  results  were  extended  by  Uchida  (Ref.  2)  to  pulsatile  flow  superimposed  on  a 
laminar  mean  fluid  motion.  Although  Uchida’s  formulation  is  valid  for  pulsations  of 
arbitrary  waveform,  his  results  consider  only  the  case  of  simple  harmonic  oscillations. 
Romie  (Ref.  3)  treated  the  velocity  and  temperature  field  in  laminar  pipe  flows  subjected  to 
simple  harmonic  pressure  pulsations  in  a  study  of  heat-transfer  mechanisms.  The  velocity 
field  in  the  des'doping  flow  near  the  entrance  of  a  tube  was  treated  by  Atabek  and  Chang 
(Ref.  10).  Barnett  (Ref.  4)  obtained  the  velocity  and  temperature  fields  in  turbulent 
pulsating  flows  in  determining  the  effect  of  simple  harmonic  longitudinal  vibrations  on  heat 
transfer  in  a  tube. 

It  is  the  purpose  of  this  study  to  obtain  a  solution  for  the  velocity  field  in  a  fully 
developed,  turbulent  pipe  flow'  subjected  to  excitation  by  longitudinal  disturbances  of 
arbitrary  waveform.  The  particular  cases  to  be  studied  will  include  (1)  simple  harmonic 
pulsations  (sinusoidal  variation  in  time)  of  the  pressure  field  or  longitudinal  vibrations  of 
the  tube,  (2)  coupled  pressure  pulsations  and  vibrations,  and  (3)  pressure  pulsations  of 
arbitrary  waveform. 

2.0  FLUID  RESPONSE  TO  LONGITUDINAL  PULSATIONS 

The  subject  of  this  study  is  the  motion  of  a  fluid  in  a  pipe  when  the  fluid  experiences 
pulsations  induced  either  by  a  time-varying  pressure  field  or  by  longitudinal  vibrations  of  the 
tube.  Two  candidate  cylindrical  reference  frames  for  formulating  the  problem  are  shown  in 
Fig.  1.  In  the  inertial  reference  frame  (r\  8' ,  z),  the  "no-slip”  condition  at  the  tube  wall 
(r  =  R)  requires  specification  of  a  nonhomogeneous,  time-dependent  velocity  at  the 
boundary  when  the  tube  is  vibrating.  If,  on  the  other  hand,  a  similarly  oriented  coordinate 
system  (r,  8,  z)  is  chosen  to  move  at  the  velocity  of  the  wall,  a  homogeneous  boundary 


W  =  -ZQu  sin  ut 


Figure  1.  Comparison  of  inertial  and  tube-fixed 
coordinate  systems. 
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condition  will  result,  and  the  problem  is  simpler  from  a  mathematical  viewpoint.  This  latter 
system,  however,  is  noninertial,  and  the  momentum  equations  must  be  modified  to  account 
for  this  fact. 

2.1  FORMULATION  IN  NONINERTIAL  FRAME 

The  continuity,  momentum,  and  energy  equations  for  the  isothermal  flow  of  an 
isotropic,  Newtonian  fluid  in  a  noninertial  coordinate  system  (Ref.  4)  are  given  below. 


Continuity: 


Dp 

m 


+  p.\ 


(D 


Momentum: 


(2) 


Energy: 


P  Di 


(3) 


where 

n  ■j  ^  a\ 

—  -  —  +  V.  —  .  A  =  -  ,  unci  $  =  dissipation  function* 

Di  J  fJx  dx 

j  j 

Comparing  these  equations  to  those  obtained  for  an  inertial  coordinate  system,  one  sees 
that  they  differ  only  by  the  term  g(dVj/dr),  which  appears  in  the  momentum  equation.  The 
vector  V,  is  the  velocity  of  the  moving  reference  relative  to  a  fixed  frame.  For  a  simple 
harmonic  oscillation  of  frequency,  qj,  and  amplitude,  Zo,  directed  along  the  Z-axis,*  the 
following  relations  may  be  established  between  the  coordinate  systems  of  Fig.  1: 

r '  =  r  ,  6'  =  0  ,  z  =  z  +  ^0  L'oS  w  t  (4) 


♦Note  that  more  complex  vibrations  may  be  treated  by  introducing  the  appropriate  reference  motion  here. 
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Thus,  if  U,  V,  W  are  the  components  of  the  coordinate  system  velocity  vector,  V„  then 

U  -  V  =-  0  and  W  =  —  Za  (u  los  tat  (5) 

The  relations  ( 1)  through  (3)  and  (5)  arc  adequate  to  describe  the  laminar  flow  of  fluid  in 
a  longitudinally  vibrating  tube.  To  obtain  a  formulation  suitable  for  turbulent  flows  as  well, 
it  is  convenient  to  decompose  any  dynamic  variable,  B  (such  as  pressure  and  velocity),  into 
the  sum  of  a  statistical  mean  value,  B  =  b,  and  a  fluctuating  component,  b'(t).  Thus: 


V.  =  v‘,v'  / 

P  -  p  -  p '  \ 


(6) 


It  should  be  noted  that  if  temporal  averages  are  utilized,  the  period,  T,  over  which  the 
average 

r 

b  =  -  f  Ril  )dl 

is  obtained  must  be  small  with  respect  to  the  period  of  the  imposed  oscillation  to  allow  the 
mean  value  to  be  time  dependent.  Introducing  Eqs.  (4)  and  (6)  into  Eqs.  (1)  through  (3), 
assuming  the  fluid  to  be  incompressible,  and  taking  the  average  of  the  resulting  equations  in 
accordance  with  Reynolds  rules  {Ref.  11)  leads  to  the  following  relations  with  respect  to  the 
noninertial  cylindrical  reference  system: 


Continuity: 

l  B  .  .  J  b  v  r)  v, 

-  —  mi)  + - —  =  0 

1  cb  r  BO  d  i 


r-Momentum 


dii  (9 u  v  f5u  v‘  dll  J  <■):> 

—  +  u  —  + - - - -*»v  —  = - 1~ 

C?t  C?r  '  B6  T  B?.  P  Br 


(7) 


(8) 
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0-Momentum: 

1  0p 

Pr  00 


0v  0v  v  dv  iiv  0v 

-  H  II  -  f - +  -  +  W  - 

01  0r  1  00  r  07. 


1  0  1 

{r*L\ 

v  1  02  V 

2  0\ 

0f_v‘ 

1 - 

^  i 

^1 

{  dr) 

'2  f2  00 2 

+  ^do  + 

0Z2_ 

a  r 
*  ,(u 


v') 


U  (v'2)  + 
r  00 


_0 

J?7 


(v  '  W  ')  -J 


7  {u'v') 

— 


(9) 


z-Momentum: 


dw 

t9w 

\  0W 

0W 

— 

f  u  — 

+  —  - 

+  w  — 

0t 

dr 

«■  00 

0z 

1  01' 
P  07. 


COS  &>L 


The  energy  equation,  Eq.  (3),  simply  stales  that  the  dissipation,  3>,  is  zero. 

These  equations  may  be  greatly  simplified  under  the  following  assumptions:  (1)  the  mean 
flow  is  rectilinear,  so  that  u  =  v  =  0;  (2)  the  mean  velocity  field  is  axisymmetric;  and  (3)  all 
correlations  (Vj'Vj1)  are  axisymmetric  and  invariant  in  the  axial  direction.  The  first  and  third 
assumptions  restrict  the  analysis  that  follows  to  the  hydrodynamically  fully  developed  region 
of  a  tube.  For  steady,  turbulent  flows  the  condition  is  usually  obtained  within  10  tube 
diameters  from  the  entrance  of  the  tube  for  Re  a:  10,000  (Ref.  12). 

For  these  assumptions,  the  continuity  equation  reduces  to  the  statement  that  w  =  w  (r, 
t).  The  momentum  equations  reduce  to: 


r-Momentum: 


0-Momentum: 


P  0r  0r 


—  (u'v')  +  -  (u'v')  =  0 

0r  r 


(ID 


(12) 
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z-Momentum: 


Qw 

Ih 


>  <3]>  ? 

- —  -I  CJ  cos  ml 

P  dr.  0 


i/  <9  /  t)w  \  ]  c) 
r  dr  \  dr  )  r  dr 


(ru'n  ') 


(13) 


The  formulation  of  the  problem  is  completed  upon  specification  of  the  boundary 
conditions.  For  a  no-slip  condition  at  the  tube  wall  it  is  required  that 


w  (R,  0  =  0  (14) 

Several  other  restrictions  must  also  be  imposed:  (1)  any  solution  must  be  symmetric  about 
r  =  0;  (2)  any  solution  must  be  finite;  and  (3)  all  transient  effects  must  vanish  as  the 
amplitude  or  frequency  of  pulsation  becomes  zero.  Since  flow  near  the  tube  boundary  is 
known  to  be  laminar,  all  correlations  must  vanish  at  the  wall.  No  initial  condition  will  be 
imposed  upon  the  problem.  Any  solution  to  be  obtained  will,  therefore,  be  quasi-steady  and 
correspond  to  conditions  existing  after  the  initial  transients  have  been  damped. 

2.2  DIMENSIONAL  ANALYSIS 

It  is  convenient  and  instructive  to  obtain  a  dimensionless  system  of  equations  equivalent 
to  those  obtained  above.  To  this  end  the  following  dimensionless  independent  variables  are 
introduced: 

r "  1  ’  - ;  :  <l5> 
The  dependent  variables  are  taken  to  be 

B(» h  £  r>  -  — — -  ;  V (jj,  ?)  --  (16) 

pi  * 2  * 

where  the  friction  velocity  is 


(Angle  brackets,  <  >,  denote  the  average  over  a  cycle  of  the  pulsation.)  The  reference 
pressure,  p*,  is  taken  to  be  the  pressure  at  the  entrance  to  the  tube.  The  correlations  are 
defined  by 


V.  V 

*  J 


(18) 
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Substitution  of  Eqs.  (15)  through  (18)  into  Eqs.  (11)  through  (13)  leads  to  the  following: 


rj-Momentum: 


3R-... 


rtrj  <h\ 


II  -  It 

II  II  V  V 

-*  -  =  0 


(19) 


0-Momentum: 


d  R  it 

u  y  _  n  \  _ 

i-  2  -  =  0 


(20) 


f-Momentum: 


Re  — 
'  * 


—  He  k  — 


\  Hcj 


I  d  (  «9V  \  R«“  d 


MS  rf-  —  [1?—  I - —  (7|  H  ) 

Uv-  V  c>t)  l  thf  /  n  07!  ,l" 


(21) 


where 


K  = 


cu  n- 

V 


is  the  dimensionless  frequency  or  vibrational  Reynolds  number, 


Re' 


u-R 

V 


is  the  friction  Reynolds  number,  and 


it 


(22) 


(23) 


(24) 


is  the  dimensionless  amplitude. 

The  amplitude  of  a  pressure  pulse  will  be  determined  later  (see  Section  2.3),  but  is  of  the 
form 


\ 

.p 


ti  I  ( i ,  =« ) 
[2  (1.  ool 

it 


(25) 
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where  fi|  {1,  oo)  is  the  amplitude  of  the  time-dependent  component  of  the  wall  pressure 
coefficient  far  from  the  entrance  of  the  tube  and  %  (1,  oo)  is  the  amplitude  of  the  steady 
flow  component. 

In  view  of  these  equations,  it  is  clear  that  the  velocity  field  in  a  fully  developed,  pulsating 
flow  is  given  by  a  relation  of  the  form 


V  V(i,,  r.  Re«,  Rev.  A,  A  ) 


(26) 


The  friction  Reynolds  number  is  related  to  the  conventional  Reynolds  number,  Re  = 
2<U>R/,V,  where  U  is  the  spatially  averaged  velocity  in  the  tube,  by 


_  Rr  _ 

8 

where  Cf  is  the  Fanning  friction  factor  (Ref.  12).  When  the  mean  velocity  <U>  =  0,  the 
friction  Reynolds  number  vanishes  and  must  be  replaced  in  the  f -momentum  equation  by 
the  harmonic  Reynolds  number  (Ref.  7): 

to  K  7. 


Re. 


A  Rr 


The  velocity,  pressure  coefficient,  and  correlations  for  no  mean  through  flow  are  thus 
nondimensionalized  with  respect  to  the  harmonic  velocity,  Uh  =  wZ0. 

2.3  SPECIFICATION  OF  THE  PRESSURE  FIELD 

Consideration  must  first  be  given  to  determining  the  form  of  the  pressure  coefficient  for 
various  types  of  pulsations.  Although  an  exact  formulation  will  not  be  obtained,  the 
functional  behavior  of  the  pressure  field  will  provide  sufficient  information  to  allow  a 
solution  of  the  momentum  equations. 

The  0-momentum  equation,  Eq.  (20),  may  be  integrated  immediately  to  obtain 


(27) 


For  the  correlation  to  vanish  at  the  tube  wall  (i?  =  1),  it  follows  that  C  =  0  and  Ruv  =  0 
everywhere  in  the  tube. 
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Direct  integration  of  the  n-momentum  equation,  Eq.  (19),  is  also  possible,  and  one 
obtains 


n  (77,  4,  r)  =  n(l,  &  r)  I  Fq  Iq,  r) 


(28) 


where 


-F(7J,  rl  =  K„  M(j r)  - 


It  -  it 

IJ II  II  V 


drj 


Note  that  while  the  pressure  varies  across  the  tube  in  a  turbulent  flow,  this  variation  is 
independent  of  the  axial  coordinate  provided  the  flow  is  fully  developed.  For  this  case, 
moreover,  it  is  reasonable  to  assume  that  at  any  given  axial  station  the  wall  pressure 
coefficient  will  exhibit  a  temporal  behavior  similar  to  the  excitation  applied  at  the  entrance 
to  the  tube.  Consequently,  for  a  simple  harmonic  pressure  oscillation, 


[° 


ft  (  I  .  0  =  C  |  C  I  ft  (I,  «C  )  -  0  1  (  ]  .  )  UO'i 


0 


Coir) 


(29) 


where  fli  (I,  <»)  is  the  mean  value  over  a  cycle  of  the  wall  pressure  coefficient  in  the  fully 
developed  region  and  fl|  (1,»)  is  the  amplitude  of  the  variation  about  the  mean.  At  f  =  0,  it 
follows  that 


C  2  ( r )  -  Q{I.O,  r)  =  n(i0,  OH-  Qjfl.  0)  cos  r 


(30) 


if  the  waveform  is  undistorted.  From  Eqs.  (29)  and  (30)  it  follows  that 


nn.  &  r)  =  (:,£«„(!, -Hi 


A  co»  r)  t  n„(],o)  (i 


i  A  cos  r) 

*  n 


(31) 


From  Eqs.  (28)  and  (31),  the  f -momentum  equation  becomes 


lie 


i9v 
v  St 


A  Ho* 

-f.0  He’  n  ( l,  «  )  ( I  l-  A  I  OH  r)  ■-  - -  (.Ob  t 

°  I’  llc“ 


1  A 

1  Off 


i-M 

drj 


11  w 


) 


(32) 
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Now  1  he  dimensionless  mass  How  through  the  tube  is 


G*  =  :I>  =  2rr  f  <V> 

. .  .  ...  '? 


pu*n2 


(33) 


If  this  value  is  to  remain  constant  over  a  cycle  of  the  pulsation,  it  follows  that  d<V>/dr  = 
0.  Accordingly,  averaging  Eq.  (32)  over  a  cycle  leads  to 


Multiplying  Eq.  (34)  by  17  and  integrating  gives 


7,  — —  -  Re1  7,  <R  >  ^  Re*  C»On  ^iC. 
1  Sr,  "*  2  o  2  i 


(35) 


Evaluating  Eq.  (35)  on  the  centerline  of  the  tube  shows  that  C4  =  0.  Furthermore,  at  the 
wall  <  Ruw  >  =  0;  thus 


C, 


n.j*  n  u,  » ) 


<kv> 

<*7 


17=1 


(36) 


From  Newton’s  law  of  viscosity  it  follows  directly  that 


d<V> 

&Tj 


So  that 


V=  1 


-HcJ 


n  ci,«) 


(37) 


The  pressure  coefficient  variation  in  the  fully  developed  region  is,  thus, 

Q(i7,£,  t)  =■  -2£U  (  A]t  cos  r)  +  J10(]T  0)^1  t-  A{)  cos  r  ^ 


(38) 
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From  Eq.  (38),  the  f -momentum  equation  becomes 


=  J(r)  +  I  £  fa-)-  —  ~  (7«  ) 

dr  *1  dr/  \  (hf  /  V  fir/  1IW' 


(39) 


where,  for  simple  harmonic  pulsations,  the  forcing  function,  .?(r),  is 

?(r>  =  2  He*  (1  T  ff  cos  T) 
with  a  generalized  amplitude  defined  by 


<?  =  A 


p 


A  He2 

v 

2  Rc*  2 


(40) 


(41) 


It  is  noted,  moreover,  that  while  Eq.  (39)  was  obtained  for  the  case  where  pressure 
pulsations  and  vibrations  of  a  common  frequency  were  the  forcing  function  for  the  transient 
fluid  motion,  the  case  of  tube  vibration  alone  (Ap  =  0)  or  pressure  pulsations  alone  (A  =  0) 
are  also  obtained.  In  fact,  any  pulsation  for  which  dQ/d?  =  J(r)  can  be  treated  by  simply 
expanding  .7(r)  in  a  complex  Fourier  series,  as  will  be  shown  in  Section  2.5. 

2.4  SIMPLE  HARMONIC  PULSATIONS 

It  has  been  shown  that  for  simple  harmonic  vibrations  of  a  tube  and/or  pressure 
pulsations,  the  velocity  field  in  axisymmetric,  fully  developed  flow  of  a  constant  property 
fluid  must  satisfy  a  differential  equation  of  the  form  ofEq.  (39).  Since  this  relation  is  linear, 
a  solution  is  sought  which  will  consist  of  a  periodic  motion  superimposed  on  the  mean  flow. 
Thus,  let 


t)  =  V0(jj)  -  V  j  (l},  T ) 


Similarly,  in  analogy  with  the  Boussinesq  exchange  hypothesis  (Ref.  11), 


It 


II W 


v  He* 


r  (iv 

<m  (l})! n 
o  □  Tf 


+  £ 


(>?) 


(42) 


(43) 


Substitution  of  Eqs.  (41)  through  (43)  into  Eq.  (39)  allows  one  to  separate  the  f -momentum 
0  equation  into 
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for  the  mean  velocity  component,  and 


(44) 


(45) 


for  the  transient  component  of  velocity.  The  boundary  conditions  corresponding  to  Eq.  (14) 
aTe  V0(  l )  =  Oand  V|(1,t)  =  0. 

2.4 J  The  Mean  Velocity  Component 

Equation  (44)  and  its  boundary  conditions  are  identical  to  the  Reynolds  equation  and 
boundary  conditions  for  fully  developed,  “steady”  turbulent  flow  in  a  tube.  Since  turbulent 
flow  solutions  rely  on  agreement  with  experimental  data  as  well  as  with  the  governing 
equations,  many  semi-empirical  solutions  to  that  problem  have  been  proposed.  Due  to  its 
inherent  simplicity,  the  solution  proposed  by  Pai  (Ref.  13)  will  be  used. 

The  velocity  variation  across  the  tube  is  given  by 


where  Vc  =  Vo(0).  This  relation  is  an  exact  solution  of  the  Reynolds  equation  if  the  eddy 
diffusiviiy  is 


m  —  s 
**  (ill  —  ] ) 


in  (*.  —  J ) 

M  (ill  —  I  ) 


r)im  - 


(47) 


The  empirical  coefficients  must  be  chosen  so  that  the  velocity  profile  accurately  represents 
the  experimental  data. 

0 
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If  one  differentiates  Eq.  (46)  and  evaluates  the  result  at  rj  =  1,  it  follows  that 


Re* 

s  -  - 

2V 

c 


(48) 


The  parameters  “s”  may  be  interpreted  as  the  ratio  of  the  wall  shear  stress  in  a  turbulent 
flow  to  that  in  a  laminar  flow  with  the  same  centerline  velocity  {Ref.  13).  For  s  =  1,  Eq.  (46) 
reduces  to  the  fully  developed  laminar  flow  profile,  and  the  eddy  diffusivity  vanishes.  To 
evaluate  “m”  it  is  only  necessary  to  require  that  the  mean  velocity  predicted  by  the  profile 
should  match  those  velocities  experimentally  observed.  Thus  (Ref.  4), 


From  these  relations,  it  is  readily  shown  that  for  any  specified  mean  velocity,  V,  both 
empirical  constants  can  be  determined  if  the  wall  shear  stress,  <  erw  > ,  and  the  centerline 
velocity,  Vc,  are  known.  Figure  2,  based  on  the  data  of  Haberstroh  and  Baldwin  (Ref.  14), 
presents  the  values  of  the  empirical  constants  for  4,000  <  Re  <  200,000.  For  Reynolds 
numbers  in  excess  of  10,000,  both  parameters  are  accurately  described  by  the  power  law 
relations  shown  in  the  figure. 

Figure  3  compares  Pai’s  relation  (Eq.  46)  to  data  obtained  by  Nikuradse  for  Re  = 
23,000.  The  agreement  is  good  over  the  entire  tube  radius  although  Pai’s  result 
underpredicts  the  experimental  data  by  about  10  percent  near  the  wall.  Figure  4  shows  the 
velocity  profile  and  eddy  diffusivity  variation  across  the  tube  for  several  Reynolds  numbers. 
It  is  noted  that  for  high  Reynolds  numbers  approaches  a  constant  value  across  the 
tube.  This  observation  will  be  of  some  importance  in  the  following  section. 

2.4.2  The  Transient  Velocity  Component 

Noting  that  cos7  =  Re  (e,r)  suggests  a  solution  to  Eq.  (45)  of  the  form 


V  |  { »j )  =  Ro 


(50) 
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Making  this  substitution  in  Eq.  (45),  one  obtains 

f  _  j  ]\c ^  r  -  -Ci  (51) 


where  primes  denote  differentiation  with  respect  to  r}.  The  boundary  condition 
corresponding  to  a  no-slip  requirement  at  the  wall  is  simple  f(l)  =  0. 

Equation  (51)  can  be  solved  by  numerical  techniques  for  any  arbitrary  eddy  viscosity 
distribution.  One  procedure  For  obtaining  such  a  solution  is  presented  in  Appendix  A.  The 
main  objection  to  such  solutions,  however,  is  a  lack  of  definitive  experimental  data  on  the 
nature  of  the  turbulent  exchange  process  in  pulsating  flows.  Bogdonoff  (Ref.  6)  concluded 


Figure  2.  Variation  of  velocity  profile  parameters  with 
Reynolds  number. 
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Figure  3.  Comparison  of  theoretical  and  experimental  velocity 
distributions  in  steady  turbulent  flow. 


Dimensionless  Radial  Location,  v 

Figure  4.  Velocity  and  eddy  viscosity  distributions  in  steady, 
fully  developed  flow  in  a  tube. 
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on  the  basis  of  hot-wire  studies  that  the  exchange  process  was  altered  but  presented  no 
quantitative  information.  Mickelson  and  Lawrence  (Ref.  15),  on  the  other  hand,  showed 
that  acoustic  excitations  affect  the  spectrum  of  turbulence  only  in  the  immediate  vicinity  of 
the  exciting  frequency.  Recent  experiments  by  Clamen  and  Menton  (Ref.  7)  show  that 
pulsating  Hows  with  no  mean  motion  begin  to  exhibit  intermittent  turbulent  behavior  when 
the  harmonic  Reynolds  number  (Reh  =  A  Rev)  exceeds  1 ,000  and  arc  fully  turbulent  for  Rctl 
>  3,000.  None  of  these  studies  permits  the  exchange  process  to  be  modeled  in  detail 
although  it  is  probable  that  turbulent  exchange  is  increased  in  the  presence  of  pulsations.  In 
the  absence  of  the  required  information,  some  insight  into  the  influence  of  pulsation- 
induced  turbulence  can  be  obtained  by  assuming  thaL  the  eddy  diffusivity  associated  with  the 
transient  velocity  field  is  spatially  invariant.  By  analogy  to  the  eddy  viscosity  distribution  in 
steady  flows  (Fig.  4),  this  assumption  should  approximate  the  physical  behavior  if  the 
vibrational  and  harmonic  Reynolds  numbers  are  large. 

For  tm[/v  =  constant,  Eq.  (51)  becomes 


where 


and 


f"  !-  -  I'  -  i  Ro*  f  ,  -ft* 
7 


lie* 


R<; 


V 


I ' 


A  particular  solution  to  Eq.  (52)  is 


f 


i* 


Hr;* 


(52) 


(53) 


(54) 


(55) 


The  homogeneous  form  of  Eq.  (52)  is  a  modified  Bessel  equation  of  order  zero.  Its  solution 
is 

fh  -  f  I  l.,0h  a‘  0  <  f'2Ko(i^  a*  rj)  (56) 
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where  a*  =  VReJ  is  commonly  termed  the  frequency  factor.  The  general  solution  to  Eq. 
(52),  then,  is  the  sum  of  Eqs.  (55)  and  (56).  For  the  boundary  condition  f(l)  =  0  and  the 
constraint  that  f(ij)  must  be  finite  for  all  values  of  tj,  0  <  t)  <  1 ,  the  general  solution  may  be 
written  as 


JL  (  ^  _!  I 

-*)  "  ( 


(57) 


It  is  convenient  to  express  this  result  in  an  alternative  form.  First,  it  is  noted  that  a 
modified  Bessel  function,  J„  (il/2x)  may  be  expressed  in  terms  of  the  Kelvin  functions  ber„x 
and  bei„x  by  (Ref.  16) 

1„  (i,!  O  =  bei'^x  +  i  bei^x  (58) 


An  equivalent  form  of  Eq.  (58)  is  the  polar  representation 

\v  ( i x)  =  (x)  exp  [  i  <?„  (x)  -  y  J 


(59) 


where  the  modulus  is 

(x)  =  ^ber2  (x)  bci2  (x)  J  ' 


and  the  phase  is 


i  beit,x 

6  (x)  =.  tan  - 

buV 


Substituting  Eq.  (57)  into  Eq.  (50)  along  with  the  definition  given  by  Eq.  (59)  gives 


ff  (  m„  ( B*  »?)  r  1 ) 

V,  (r),  t)  =  -  <  sin  r - - —  sin  r  -  0  (  a*  jf)  -  6  (  a*)  > 

1  lie  )  M  (a*)  L  °  ° 

v  ",  n  i 


(60) 


Detailed  discussion  of  this  result  will  be  deferred  until  Section  3.0,  but  two  effects  are  readily 
apparent  from  the  equation.  First,  the  magnitude  of  the  transient  velocity  is  directly 
proportional  to  an  attenuated  value  of  the  generalized  amplitude.  For  a  fixed  value  of(S, 
high  frequencies  tend  to  diminish  the  effect  of  pulsations.  Secondly,  the  fluid  velocity  is  out 
of  phase  with  the  excitation,  S  cos  r,  by  an  amount  which  varies  with  location  in  the  tube. 
On  the  centerline  of  the  tube,  it  is  readily  apparent  that  the  phase  lag  approaches  90  deg. 
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2.5  PULSATIONS  OF  ARBITRARY  WAVEFORM 

In  many  cases  of  engineering  significance,  (he  forcing  function,  Jfr),  of  Eq.  (39)  is  a 
periodic  function  of  greater  complexity  than  the  simple  harmonic  pulsation  discussed  above. 
If,  however,  the  function  is  bounded  and  has  a  finite  number  of  discontinuities  in  the 
interval  P  =  2tt/w,  it  is  known  that  a  Fourier  trigonometric  series  will  converge  to  the  value 
of  the  forcing  function  in  the  interval.  Consequently,  many  practical  waveforms  can  be 
represented  by  a  Fourier  series  in  the  various  harmonics  or  the  basic  frequency,  w. 
(Turbulent  fluctuations,  however,  are  aperiodic  and  require  a  Fourier  integral  approach.) 

Considering  the  approach  of  the  preceding  section,  it  is  desirable  to  express  the  forcing 
function  by  the  complex  scries 


3<t)  =  ft. 


I 

n-  | 


ft. 


(61) 


where 


ft 


n 


ft  ( r') 


dr' 


(62) 


In  Appendix  B  it  is  shown  that  this  expansion  is  a  special  case  of  the  more  general  complex 
Fourier  expansion.  It  applies  for  any  J(t)  which  is  real  and  satisfies  the  Dirichlel  conditions. 
The  limits  of  integration  in  Eq.  (62)  extend  from  an  arbitrary  point  C  of  a  cycle  to  the  initial 
point  of  the  succeeding  cycle  while  the  summation  indices  are  harmonics  of  the  basic 
frequency  of  the  oscillation,  w  =  27r/P. 

A  form  of  the  velocity  field  is  sought  which  is  consistent  with  Eq.  (61).  Accordingly,  let 


V  (tj,  r)  -  V0  +• 


£  Vn(7)el 

n  -- 1 


Similarly,  the  correlation  function  may  be  written  as 


V  fll 


<1 V 


dr) 


X ,,)  ^ 


X 


(63) 


(64) 
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Substitution  of  Eqs.  (61),  (63),  and  (64)  into  the  axial  momentum  Equation,  Eq.  (39), 
leads  to  a  separation  of  variables  in  which  the  following  relations  must  be  satisfied: 


Equation  (65)  governs  the  mean  fluid  motion  and  is  identical  to  Eq.  (44)  for  Re*  = 
Jf„/2.  Since  the  boundary  conditions  to  be  satisfied  are  V0(l)  =  0  and  Vn(l)  =  0,  the 
solutions  proposed  in  Section  2.4.1  for  the  mean  motion  arc  applicable  here  as  well.  For  Eq. 
(66)  to  be  valid  at  all  times,  it  follows  that  each  harmonic  must  obey  a  relation  of  the  form 


I  i  r  (  \  rIV 

'  —  rit  \  h  -  )  — -i  n  Rc  V 

V  d7l  I  l'  /  <17,  J 


—2  K_ 


or 


(■ 


V 


i  n  Rey 


-2  Hn 


(67) 


The  similarity  of  Eqs.  (67)  and  (5 1)  is  readily  apparent,  as  is  the  observation  that  again  a  lack 
of  information  exists  on  the  eddy  diffusivity  associated  with  the  transient  fluid  motion.  As 
before,  the  heuristic  assumption  will  be  made  that  the  eddy  viscosity,  em],is  spatially 
invariant.  The  solution  to  Eq.  (67)  is  then 


2  i  H  \  M  (  a*  7i  I 

n  }  n  ft  1 

■  ■  ■  ■  -  ■  \  — . . — - 

11  /  M  (a*  ) 

i  o  n 


<■ a ;  7  >  -  e0 1 a 


(68) 
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where 


(69) 


Inlroducing  Eqs.  (68)  and  (46)  into  (63)  completely  specifies  the  velocity  field  in  the  tube 
at  any  time.  For  computational  purposes  it  is  expedient  to  compute  the  elements  of  Eq.  (63) 
individually  and  sum  the  results,  Additional  savings  in  computational  time  are  obtained  if 
the  various  coefficients  can  be  related  to  the  coefficients  of  a  Fourier  trigonometric  series 
(see  Appendix  B)  so  that  a  Fast  Fourier  Transform  (FFT)  algorithm  (Ref.  17)  may  be  used. 
With  this  in  mind,  note  that  Eq.  (62)  may  be  rewritten  as 


2Kn 


(70) 


where 


A 


r. 


If(r')  cos  nr'  dr' 


=  0,  ! ,  2 


(71) 


and 
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From  Eq.  (70)  and  DeMoivre’s  relation 

e"1  =  cos  x  -r  i  sin  x 

it  can  be  shown  by  a  straightforward  substitution  that  the  real  part  of  the  nth  harmonic  of  the 
velocity  may  be  expressed  as 
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where 


;\n  (r)  =  cos  nr  -  An  sin  nr 


(74) 
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An  equivalent  result  to  Eq.  (73)  is 
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As  in  the  case  of  simple  harmonic  oscillations  the  amplitude  of  the  transient  component  is 
attenuated  at  high  frequencies,  and  phase  lags  exist  that  depend  on  radial  location  and 
frequency,  A  more  detailed  discussion  of  these  results  and,  in  particular,  the  effect  of 
waveform  on  the  solution,  is  presented  in  Section  4,0. 

3.0  THE  VELOCITY  FIELD  FOR  SIMPLE  HARMONIC  PULSATIONS 

The  response  of  the  fluid  to  a  pulsation  of  the  form  described  by  Eq.  (40)  is  simply  the 
sum  of  Eqs.  (46)  and  (60).  Although  this  sum  is  a  closed  form  analytic  solution  of  the 
transient  Reynolds  equation  for  fully  developed  flow,  the  complexity  of  the  equations  and 
the  number  of  auxiliary  relations  to  be  solved  make  its  evaluation  practical  only  by  use  of  a 
computer.  A  program  identified  as  PLJLSAT  has  been  written  for  an  IBM  370/165  digital 
computer  to  evaluate  the  velocity  field  in  fully  developed,  laminar  or  turbulent,  steady  or 
pulsating  flows  of  simple  harmonic  waveform. 

3.1  THEORETICAL  RESULTS 

The  radial  distribution  of  velocities  in  a  flow  with  a  mean  Reynolds  number  of  10,000  is 
shown  in  Fig.  5  for  several  limes  during  a  cycle  of  pulsation.  At  t  =  0,  the  velocity  profile  is 
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seen  to  correspond  to  steady  flow  in  a  tube  at  this  Reynolds  number  except  for  a  thin  layer  of 
fluid  near  the  tube  boundary  where  the  velocity  increases  rapidly  to  a  value  in  excess  of  that 
observed  on  the  centerline  before  decaying  to  zero  at  the  wall.  It  should  be  recalled  that,  due 
to  the  quasi-steady  nature  of  the  solution,  t  —  0  corresponds  to  the  start  of  a  cycle  but  not  to 
the  initial  time  at  which  pulsations  are  induced.  With  increasing  values  of  r  the  magnitude  of 
the  velocity  increases  until  the  maximum  value  is  obtained  at  t  -  90 deg.  The  magnitude  and 
extent  oT  the  velocity  increases  near  the  wall  are  also  seen  to  vary  with  time  during  a  cycle. 
After  the  maximum  velocity  is  obtained,  the  velocity  is  seen  to  decrease  until  the  minimum 
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Figure  5.  Temporal  variation  of  local  velocity 
for  harmonic  forcing  function. 
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values  are  obtained  at  r  =  270  deg.  It  is  apparent  and  will  be  seen  in  later  results  that  during 
a  cycle  there  are  times  at  which  the  core  velocities  are  positive  while  near  the  wall  the  velocity 
is  negative. 

Since  the  forcing  function  is  a  maximum  at  r  =  0  and  a  minimum  at  t  =  180  deg,  the 
centerline  velocity  lags  the  forcing  function  by  90  deg.  This  behavior  was  anticipated  on  the 
basis  of  Eq.  {60),  which  shows  that  the  phase  lag  depends  on  both  the  frequency  parameter, 
a*  =  V Rev,  and  the  radial  location,  ij.  In  Fig.  6,  the  phase  lag  is  plotted  as  a  function  of  j? 
for  two  vibrational  Reynolds  numbers.  The  fluid  phase  behavior  is  of  significance  in  any 
application  in  which  a  secondary  response  to  the  forcing  function  (such  as  particle  dynamic 
effects,  flow-induced  vibrations,  etc.)  are  of  importance.  A  striking  example  is  that  some 
techniques  proposed  for  sizing  particles  contained  in  a  flow  (Ref.  18)  relate  the  particle  drag 
to  pressure  pulsations  imposed  upon  the  flow.  Failure  to  properly  account  for  the  fluid 
behavior  could  lead  to  serious  errors  in  such  determinations. 


Figure  6.  Phase  lag  of  the  velocity  for  a 
simple  harmonic  pulsation. 

The  effect  of  the  mean  throughflow  Reynolds  number  on  the  velocity  observed  at 
various  times  during  a  cycle  in  which  pulsations  arc  vibration-induced  is  illustrated  in  Fig.  7. 
As  would  be  expected,  the  gross  effect  of  increasing  Reynolds  number  is  to  diminish  the 
overall  influence  of  the  pulsation.  Indeed,  for  a  fixed  value  of  A  and  ReCl  the  transient 
component  of  velocity,  Eq.  (60),  is  independent  of  Reynolds  number  provided  the  turbulent 
structure  is  unaltered  or  varies  independently  of  Reynolds  number. 
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Radial  Location,  17 

Figure  7.  Temporal  variation  of  velocity  for 
several  Reynolds  numbers. 


Specifying  the  amplitude  of  the  vibration  at  the  lube  boundary.  A,  is  not,  however, 
equivalent  to  specifying  a  constant  generalized  amplitude,  as  may  be  seen  from  Eq.  (41). 
Since  the  transient  velocity  component  profiles  normalized  with  respect  to  the  centerline 
values  Vi(0,t)  are  identical  for  all  Reynolds  numbers  and  vibrational  Reynolds  numbers,  the 
effect  of  amplitude  may  be  studied  by  simply  considering  the  magnitude  of  the  centerline 
velocity.  This  evaluation  is  facilitated  if  the  radial  distribution  fuention  of  Eq.  (50)  is 
written  as 
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f(r?)  =  fr(>j  )  +  i  r  j  (77 ) 
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where,  from  Eq.  (57), 
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and 
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The  transient  velocity  component  is  then,  simply, 


V  j  (t),  *)  =  ff  cos  t  —  f1  sin  r 


(81) 


Various  authors  (Refs.  3  and  6)  define  the  strength  of  pulsation  as  the  ratio  of  the 
maximum  velocity  amplitude  on  the  centerline  of  the  tube  to  the  mean  through  flow  velocity. 
In  the  present  notation,  this  may  be  written  as 

V.  (0,  r*) 

c '  1  _ 


where  t*  is  the  time  at  which  the  maximum  value  of  the  transient  occurs.  From  Eq.  (81), 
V|  Cj?,t)  attains  a  maximum  value  for 
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The  corresponding  maximum  centerline  value  can  be  obtained  as  (Ref.  4) 
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For  all  but  extremely  low  frequencies,  the  modulus  M0(a*)  is  extremely  large,  and  the 
bracketed  term  in  Eq.  (83)  approaches  unity.  From  the  definition  of  the  friction  Reynolds 
number  and  the  mean  throughflow  Reynolds  number  the  strength  of  pulsation  is  then 
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From  Eq.  (41)  it  is  readily  apparent  that  a  constant  pulsation  strength  is  obtained  in  a 
vibrating  tube  (Ap  =  0)  for 
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Similarly,  the  strength  of  pulsation  for  pressure-induced  oscillations  (A  =  0)  is 
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Since  the  friction  Reynolds  number  is  related  to  the  mean  Reynolds  number  through  the 
Fanning  friction  factor  (see  Section  2.2),  one  of  the  parameters  can  be  eliminated  from  these 
equations.  For  example,  for  turbulent  flow  in  a  pipe  (Ref.  12), 

C,  -  0.0-16  Rc  ,'5  (Rf  <  I05) 
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for  vibrations,  and 
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for  pressure  pulsations.  It  follows,  therefore,  that  the  effect  of-  pulsations  is  dimished  at 
large  Reynolds  numbers  for  either  pressure-induced  or  vibration-induced  oscillations 
although  this  trend  is  much  more  pronounced  in  the  latter  case.  On  the  other  hand,  high 
frequencies  increase  the  pulsation  strength  for  vibrations  but  attenuate  the  effect  of  pressure 
pulsations. 

The  effect  of  frequency  on  the  velocity  profile  in  the  vicinity  of  the  tube  wall  is  shown  for 
several  vibrational  Reynolds  numbers  in  Fig.  8.  Since  these  profiles,  as  well  as  the  centerline 


Figure  8.  Influence  of  vibrational  Reynolds  number 
on  transient  velocity  component. 
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velocity  (Fig.  9)  are,  ostensibly,  independent  of  the  mean  flow  characteristics,  the  trends  are 
applicable  for  all  Reynolds  numbers.  It  is  seen  that  as  the  frequency  is  increased  the  local 
maximum  “moves”  nearer  the  wall.  This  phenomenon  has  been  termed  the  “annular” 


Figure  9.  Temporal  and  amplitude  dependence  of  transient 
centerline  velocity. 


effect  and  has  long  been  of  interest  to  fluid  dynamicists  (Refs.  1  and  9).  The  actual  location 
and  magnitude  of  the  maximum  during  a  cycle  is,  however,  a  function  of  time,  as  could  be 
inferred  from  Figs.  5  and  6.  To  illustrate  the  frequency  dependence  of  the  maximum  velocity 
location,  it  is  useful  to  consider  the  root-mean-square  (rms)  value  of  the  transient  velocity. 
From  Eq.  (81),  this  is  simply 
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The  variation  of  the  rms  velocity  is  shown  in  Fig.  10  for  a  range  of  vibrational  Reynolds 
numbers  at  a  mean  flow  of  Re  =  104.  It  is  seen  that,  for  this  vibration-induced  example,  the 
rms  velocities  increase  by  nearly  three  orders  of  magnitude  when  the  frequency  is  varied  by  a 
factor  of  20.  (Note  that  for  water  in  a  1-in. -radius  rube,  the  corresponding  circular 


Radial  Location,  rj 

Figure  10.  Effect  of  vibrational  Reynolds  number  on 
RMS  velocity  component. 


frequencies  are  0. 1  Hz  at  Rev  =  100  and  20  Hz  at  Rev  =  20,000.  For  air,  Rev  =  100  implies 
that  f  =  0.4  Hz.)  The  annular  effect  is  clearly  seen,  and  the  maximum  velocity  location,  tj*, 
is  very  near  the  wall  at  large  values  of  Rev  (Fig.  11).  It  can  be  shown,  moreover  (Ref.  19), that 
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(90) 


This  result  is  of  importance  in  calculations  of  pulsating  flow  characteristics  since  (1) 
important  details  of  the  flow  may  be  lost  if  the  computational  grid  is  not  varied  as  the  wall  is 
approached,  and  (2)  alterations  of  the  turbulent  exchange  process  may  be  related  to  the 
location  of  these  maxima  (Ref.  4). 


Vibrational  Reynolds  Number,  Rey 

Figure  11.  Variation  of  transient  component  maximum 
value  location  with  frequency. 

Both  the  analysis  and  the  PULSAT  program  allow'  the  determination  of  the  velocity 
profile  for  the  case  where  a  constant  eddy  viscosity  is  associated  with  the  velocity  transient. 
While  decidedly  heuristic,  the  assumption  leads  to  some  interesting  results,  as  may  be  seen  in 
Fig.  12.  The  velocity  profiles  showm  during  the  first  half-cycle  of  a  pulsation  have  been 
computed  on  the  basis  that  (1)  only  molecular  diffusion  occurs  (em|  =  0)  and  (2)  the 
turbulent  diffusivity  is  equal  to  the  spatial  average  of  the  eddy  viscosity  associated  with  the 
mean  How,  €m,  =  A  decided  effect  of  em }  =£  Ois  observed  at  all  times.  In  general,  since 
the  governing  parameter  for  the  altered  case  is 
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the  solution  behaves  as  though  it  corresponds  to  an  e1T1 ,  =  0  flow  at  lower  frequencies.  The 
phase  shift,  moreover,  depends  on  ReJ1'2,  and  thus  events  occur  at  different  times  in  the 
cycle  than  would  be  expected  if  alteration  of  the  diffusion  mechanisms  were  not  present. 


Figure  12.  IMear-wall  velocity  profiles  for  turbulent 
exchange  modification. 

3-2  COMPARISONS  WITH  EXPERIMENT 

Data  on  the  detailed  variation  of  the  velocity  field  in  pulsating  flows  are  extremely 
limited  although  considerable  experimental  information  on  the  effect  of  pulsations  on  heat 
transfer  exists  (Ref.  8).  Richardson  (Refs.  1  and  20)  measured  the  rms  velocity  field  in  a 
circular  tube  in  the  absence  of  a  mean  flow  component  and  identified  the  annular  effect 
discussed  earlier.  Additional  hot-wire  studies  were  made  by  Bogdonoff  (Ref.  6)  and 
Mohajery  (ReT.  21)  for  air  flows  at  mean  Reynolds  numbers  in  the  range  from  50,000  to 
100,000.  Both  investigations  showed  that  pressure  pulsations  altered  the  eddy  viscosity 
distribution  in  the  flow,  but  whereas  Mohajery  concluded  the  effect  was  not  appreciable  at 
any  radial  location  in  the  lube,  Bogdonoff  observed  significant  increases  in  eddy  viscosity 
near  the  tube  wall.  In  the  latter  study,  however,  the  author  observed  that  this  was  a  tentative 
result  since  it  was  difficult  to  compensate  for  probe  blockage  and  compressibility  effects. 
Recently,  Clamen  and  Menton  (Ref.  7)  used  a  hydrogen  bubble  technique  to  study  water 
flows  in  a  vibrating  tube  from  Re  =  0  to  Re  =  2,900  for  a  range  of  amplitudes  and 
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frequencies.  Since  these  data  conform  better  to  the  incompressible  assumption  of  the  present 
analysis,  the  comparisons  that  follow  will  be  based  on  their  data. 

Figure  13a  compares  the  velocity  distribution  measured  in  a  purely  oscillatory  flow  (Re 
=  0)  to  the  present  theory  for  several  times  during  a  cycle.  Clearly  the  theory  and  experiment 
are  in  good  agreement  at  the  times  shown  with  the  exception  of  the  small  amplitude 
oscillaiion  for  r  =  30  deg.  Since  the  hydrogen  bubble  technique  requires  synchronization  of 
a  Hash  photograph  with  a  particular  cycle  time,  it  is  possible  that  the  discrepancy  is  due  to  a 
timing  error.  (For  w  =  1 .2  Hz,  a  0.1-sec  error  would  readily  account  for  the  difference.)  In 
Fig.  13b,  theory  and  experiment  are  again  compared  for  a  mean  flow  Reynolds  number  in 
the  laminar  How  regime  (Re  =  1,535).  For  both  cases  the  agreement  is  excellent. 

A  final  comparison  is  made  in  Fig.  14  for  a  turbulent  mean  now  at  Re  =  2,900.  At  r  = 
0,  the  present  theory  represents  the  now  accurately  and  represents  a  decided  improvement 
over  a  laminar  now  theory  used  in  Ref.  7.  For  later  times  in  the  cycle,  however,  the  turbulent 
theory  grossly  overpredicts  the  measured  velocity.  While  this  may  be  due  to  a  deficiency  of 
the  theory,  it  is  more  likely  due  to  one  or  more  experimental  factors.  Clamen  and  Menton 


Figure  13.  Comparison  of  experiment  and  theory  in 
laminar  pulsating  flow. 
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note  that  at  turbulent  Reynolds  numbers  the  flow  was  “highly  disturbed”  so  that  the 
photographic  representation  of  the  velocity  profile  was  smeared.  Under  those  conditions  the 
hydrogen  bubble  technique  shows  a  range  of  velocities  so  that  the  values  obtained  are  highly 
dependent  upon  the  fashion  in  which  the  data  were  reduced.  The  possibility  or  timing  errors 
also  exists,  as  was  noted  in  the  discussion  of  Fig.  13a.  For  w  =  1.12  rad/sec,  the 
displacement  in  centerline  velocity  from  the  theory  could  be  accounted  for  by  a  time  delay  of 
0.4  and  0.2  sec  for  t  =  90  and  210  deg,  respectively.  While  such  shifts  seem  excessive,  some 
credence  is  lent  to  the  possibility  by  comparing  the  theoretical  results  for  t  =  60  deg  to  the 
data  for  t  —  90  deg.  The  agreement  is  excellent  for  ty  <  0.80. 

Re  =  2,  900 
A  =  5. 80 
Rev  =  712. 8 

- Theory,  em^  »  0 


Figure  14.  Comparison  of  experiment  and  theory 
in  turbulent  pulsating  flow. 
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Since  the  model  developed  in  Section  2.4  allows  the  investigation  of  a  uniform  alteration 
of  the  turbulent  structure  of  the  flow,  such  a  comparison  is  made  in  Fig.  14.  Clearly  the 
assumption  that  emi  =  0  is  adequate  at  t  =  0.  At  t  =  90  deg  (or  60  deg),  however,  emi  =  0 
leads  to  an  overprediction  of  the  data  obtained  near  the  wall.  By  letting  em,  =  3e^,  the 
velocity  profile  is  damped  relative  to  the  undisturbed  (em|  =  0)  profile  for  the  wall  region. 
The  results  are  in  error  at  intermediate  locations  0.4  <  ij  <  0.85,  however,  so  that  an 
assumption  of  =  constant  fails  to  predict  the  velocity  distribution  across  the  tube.  From 
these  Limited  data  it  appears  that  an  accurate  model  of  the  alteration  of  turbulence  structure 
must  take  into  account  the  radial  variation  of  the  exchange  coefficient. 

4.0  THE  VELOCITY  FIELD  FOR  ARBITRARY  WAVEFORMS 

The  analysis  of  Section  2.5  shows  that  the  response  of  the  fluid  velocity  to  a  forcing 
function  of  arbitrary  waveform  can  be  obtained  by  expanding  the  function  in  a  complex 
Fourier  series.  Each  harmonic  of  the  transient  velocity  can  then  be  obtained  in  terms  of  the 
various  harmonics  of  the  pressure  disturbance.  The  amplitude  of  the  velocity  harmonics  is 
determined  by  the  Fourier  coefficients  required  to  describe  the  forcing  function.  A  computer 
program  (ARBPLLSAT)  for  evaluating  the  velocity  response  has  been  written  for 
performing  such  calculations.  The  results  of  computations  performed  via  ARBPULSAT  for 
several  waveforms  are  described  in  the  remainder  of  this  section.  No  attempt  has  been  made 
to  comprehensively  study  the  effect  of  waveform  on  the  velocity  distribution  in  a  tube,  but 
such  studies  are  within  the  capabilities  of  the  program. 

4.1  EFFECT  OF  WAVEFORM 

To  study  the  possible  influence  of  waveform  upon  the  velocity  distribution  in  pulsating 
flows,  four  different  types  of  symmetric  waves  were  considered.  In  each  case,  the  initial  and 
final  third  of  a  cycle  were  considered  to  be  at  a  constant  pressure  while  the  pressure  varied  in 
the  central  portion  of  a  cycle  to  produce  (1)  an  instantaneous  pressure  change  followed  by  a 
later  instantaneous  decay,  (2)  a  linear  rise  followed  by  later  linear  decay,  (3)  a  linear  rise  to  a 
maximum  value  followed  immediately  by  a  linear  decrease,  and  (4)  a  half-cycle  sinusoidal 
variation  from  the  initial  pressure.  These  pressure  variations  will  be  called,  respectively,  (I)  a 
square  wave,  (2)  a  trapezoidal  w'ave,  (3)  a  triangular  wave,  and  (4)  a  sinusoidal  wave  in  the 
following  discussion. 

Figure  15  shows  the  four  waves  which  have  been  chosen  so  that  in  every  case  the  mean 
pressure  coefficient  over  a  cycle  corresponds  to  a  time-averaged  Reynolds  number  of  104. 
Three  amplitudes  are  shown  in  each  figure,  corresponding  to  peak-to-undisturbed  pressure 
amplitudes  of  fl[/Q0  =  100,  10,  and  1.5.  Also  shown  in  the  figures  are  the  Fourier  series 
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representations  of  the  maximum  amplitude  wave  for  2,048  (2n)  sine  and  cosine  transforms, 
Eqs.  (71)  and  (72).  In  all  cases  except  the  square  wave  the  greatest  departure  of  the  computed 
value,  flF,  from  the  input  function,  fy,  is  |Q|  -  £V|  '<  1.0.  Since  the  minimum  value  of 
fl[  »  25,  the  error  in  representing  the  function  at  any  point  of  the  cycle  is  less  than  4  percent 
and  is  typically  better  than  1  percent.  Even  better  accuracy  was  obtained  for  lower  values  of 
the  amplitude.  The  square  wave  is  discontinuous  at  r  =  120  and  240  deg  and  converges  to 
1/2  (J2i  +  fl0).  Since  the  computed  limit  differs  from  this  value  by  approximately  1/3  it  is 
clear  that  more  coefficients  are  necessary  to  accurately  describe  the  square  wave  in  the 
immediate  vicinity  of  the  discontinuity.  Since  this  is  a  localized  error,  it  was  not  deemed  of 
sufficient  importance  for  the  present  study  to  warrant  the  additional  computation  time. 
(Since  an  FFT  routine  computes  2n  coefficients,  the  next  possible  choice  was  to  double  the 
number  of  coefficients  used.) 

N 

The  centerline  transient  velocities,  V|(0,r)  =  E  Vn(0,T),  corresponding  to  the  various 

’]  a  I 

waveforms  are  shown  in  Fig.  16.  It  is  seen  that  regardless  of  waveform  the  centerline  velocity 
histories  are  strikingly  similar  in  both  shape  and  amplitude.  The  observed  “sawtooth” 
history  was  “somewhat  unexpected  and  immediately  raised  a  question  as  to  whether  the 
velocity  field  converges  as  rapidly  as  the  series  representation  of  the  pressure  transient.  To 
explore  this  question  it  is  useful  to  rewrite  the  forcing  function  as 


ff(-)  =  T0  +  ?,  in 


(91) 


From  Eqs.  (61)  and  (70),  then,  one  can  write 


or 


.f  I  f  r  1  (A  n  c<»*  nr  -  Hn  sin  nr) 

n^l 


C„COa  [flr  ~  dn' 


(92) 


where 


Cn  -  \IA n  "  B„2  I 


(93) 
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Figure  16.  Centerline  velocity  response  to  symmetric  waves. 
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On  the  centerline  of  the  tube  a  similar  result  can  be  obtained  if  it  is  noted  that 
M0(an*)  >  >  >  l  for  all  but  extremely  low  frequencies.  Thus  for  Rcv  >  1,  Eq.  (77)  reduces 
to 

N  c 

vV°’r)=  2  (94) 

where  Cn  and  are  again  given  by  Eq.  (93).  Clearly  this  result  converges  more  rapidly  than 
,‘fi(r)  so  that  for  any  identical  number  of  terms  the  series  For  V|(0,r)  will  provide  a  better 
representation  of  the  result  than  does  that  for  F,(t).  The  expected  rapid  convergence  of 
Vj (17.T)  is  substantiated  by  Fig.  17  where  the  transient  velocity  corresponding  to  the  partial 
sums,  E  Vn(fj,r)  is  compared  to  the  result  obtained  for  N  =  2".  Convergence  was  least 
rapid  at  r  =  120  deg,  which  is  where  the  pressure  pulse  representation  was  least  accurate. 
Even  in  this  “worst”  case,  the  error  for  as  few  as  50  terms  of  the  series  was  less  than  0.65 
percent  and  was  virtually  unaffected  by  radial  location.  It  is  concluded,  therefore,  that  the 
velocity  response  observed  on  the  centerline  of  the  tube  is  an  analytically  valid  result  and  is 
not  due  to  any  computational  deficiency  of  the  scries  representation. 


Figure  17.  Convergence  of  transient  velocity  series 
for  square  wave  excitation. 
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Some  insight  into  the  nature  of  the  behavior  or  the  centerline  velocity  is  gained  if  the 
following  integral  representation  is  considered.  Left  f(r)  be  an  arbitrary  function 
representing  the  variation  of  the  pressure  around  the  mean  and  consider  the  integral 
transform 


»<')  -  go  ’“//fir)  dr 

If  f(t)  =  a  cos  t,  a  simple  harmonic  forcing  function,  Eq.  (95)  becomes 


(95) 


g  (  - )  =  an  uos  r 


(96) 


From  Eq.  (60),  the  fluid  response  on  the  centerline  of  the  tube  is  then 


V  ,  < 0.  r) 


a  mii  r 


lib 


V 


(97) 


so  that  the  two  results  are  identical  if  a  =  Rev  ~ 1  and  ga  =  0.  For  the  square  wave.  Fig.  15, 
f(r)  =  l/3(fi0  -  ft|)  for  the  initial  and  final  portions  of  the  cycle  while  for  2ir/3  <  t  <  4n/3, 
f(r)  is  equal  to  2/3  (ft,  -  ftG).  The  transform  thus  is  of  the  form 


g(r)  =  g„  -i  a  K  r  (98) 

so  that  the  centerline  velocity  varies  linearly  in  each  region  of  the  pulse.  The  slope  aK  is, 
however,  negative  in  the  initial  and  final  stages  since  Sl0  is  less  than  Q|  and  is  positive  in  the 
central  region.  The  initial  and  final  regions  of  the  other  waveforms  also  lead  to  a  linearly 
decreasing  function,  but  the  response  to  the  trapezoidal  wave  is  parabolic  during  the  linear 
rise  and  decay  at  the  beginning  and  end  of  the  pressure  pulse.  The  triangular  and  sinusoidal 
waves,  moreover,  vary  in  a  parabolic  and  cosinusoidal  fashion,  respectively,  throughout  the 
central  region  of  the  pressure  pulse.  The  centerline  transient  velocity  in  a  pulsating  flow  is 
accordingly  given  by  a  simple  integral  transform  of  the  excitation  for  all  waveforms.  This 
behavior  is  observed  at  most  radial  locations  in  the  tube  so  that  molecular  or  turbulent 
diffusion  is  significant  only  in  the  vicinity  of  the  tube  wall.  This  is  borne  out  by  Fig.  18, 
where  the  fluid  response  to  a  square  wave  is  shown  at  several  radial  locations  in  the  tube.  For 
17  <  0.90,  the  velocity  history  is  nearly  independent  of  location,  but  near  the  wall  significant 
departures  both  in  the  magnitude  and  waveform  are  observed.  It  is  concluded,  therefore, 
that  viscous  interactions  become  important  near  the  boundary. 

It  is  clear  also  that  phase  shifts  are  characteristic  of  the  fluid  response  for  the  various 
waveforms  studied.  If  one  defines  the  phase  shift  as  the  time  delay  between  the  attainment 
of  maximum  pressure  and  velocity,  a  comparison  of  Figs.  15  and  16  shows  that  the 
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centerline  phase  lag  varies  from  40  deg  for  the  sinusoidal  variation  to  120  deg  for  the  square 
wave.  It  is  apparent  from  Figure  18,  moreover,  that  the  phase  lag  is  a  function  of  radial 
location,  as  was  shown  for  simple  harmonic  oscillations. 


Figure  18.  Comparison  of  velocity  response  to  square  wave 
at  several  radial  locations. 

The  radial  distribution  of  the  transient  velocity  component  tor  square  and  triangular 
waves  is  shown  in  Fig.  19.  The  annular  flow  phenomenon  noted  for  harmonic  oscillations  is 
retained  in  both  cases,  and  the  magnitude  of  the  maximum  velocity  is  relatively  insensitive  to 
waveform  for  all  excitations  studied.  A  slight  variation  in  the  location  of  the  maximum 
velocity  point  docs  result,  however,  for  different  forcing  functions.  Figure  20  illustrates  the 
elTect  of  vibrational  Reynolds  number  on  the  velocity  distribution  obtained  in  response  to  a 
trapezoidal  waveform.  Clearly,  increasing  frequencies  move  the  maximum  velocity  point 
nearer  the  wall  and  also  decrease  the  amplitude  of  the  velocity  transient.  Both  effects  were 
also  observed  in  Section  3.0  for  a  simple  harmonic  oscillation.  The  lowest  frequency  case 
(Rcv  =  800)  compares  the  square  and  triangulai  wave  results  to  those  obtained  for  the 
trapezoid  and  emphasizes  that  lor  comparable  amplitudes  (S2]/S20)  the  results  are  very 
insensitive  to  waveform. 
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Figure  19.  Velocity  distribution  for  various  times 
during  a  cycle. 

Il  was  noted  in  Section  3.1  that  a  convenient  measure  of  the  gross  effect  of  a  pulsation 
was  afforded  by  the  strength  of  pulsation,  Eq.  (82).  The  complex  waveforms  used  in  this 
section  do  not  lend  themselves  to  a  simple  evaluation  of  pulsation  strength  since 
each  term  of  Eq.  (94)  achieves  its  maximum  independently  of  the  others  when 
T*  =  1/n  ([(2m  -l)/2]  +  <)>„} .  For  the  higher  harmonics,  then,  r*  is  approximately  equal  to 

zero,  but  the  phase  lag  of  the  lower  frequencies  must  be  evaluated  individually  to  determine 
the  time  at  which  Vj(0,r)  is  a  maximum.  This  result  is  obviously  dependent  upon  the 
waveform  considered.  For  the  waves  considered  in  this  section,,  however,  the  centerline 
velocities  achieve  a  maximum  at  nearly  the  same  time  in  the  cycle  and  have  nearly  the  same 
amplitude.  Since,  as  can  be  seen  in  Eq.  (94),  only  the  lower  harmonics  contribute 
significantly  to  V  ](0,t),  it  seems  apparent  that  the  fluid  response  is  dictated  only  by  the  lower 
frequency  components  of  the  series. 


45 


A  E  DC-TR  -80-3 1 


Trapezoidal  Wave 


o  Square  Wave 
a  Triangular  Wave 


Figure  20.  Effect  of  frequency  on  amplitude  and  radial 
distribution  of  transient  velocity  component. 

4.2  Fluid  Response  to  a  Fluttering  Valve 

A  final  illustration  of  the  utility  of  the  computer  program  developed  in  this  study  is 
afforded  by  considering  a  basically  steady  flow  that  is  perturbed  at  regular  intervals.  Such  a 
situation  could  conceivably  be  caused  by  a  fluttering  valve  or  a  faulty  pressure  regulator 
located  in  a  flow  system.  Three  illustrative  cases  are  shown  in  Fig.  21  where  a  0.1-sec 
pressure  excursion  occurs  asymmetrically  one,  two,  or  three  times  per  second.  As  in  the 
preceding  section,  the  mean  flow  Reynolds  number  is  maintained  at  Re  =  104;  thus  the 
relative  amplitude  of  the  pulse  decreases  as  the  repetition  rate  rises. 
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Figure  21.  Asymmetric  wave  forms. 


The  fluid  response  on  ihe  cemerline  of  the  tube  is  shown  in  Fig.  22.  For  a  single 
repetition  the  velocity  history  is  similar  to  the  square  wave  siudied  earlier,  but  as  the 
repetitions  are  increased,  the  velocity  history  increases  in  complexity.  For  each  additional 
repetition  an  additional  local  maximum  (corresponding  to  the  trailing  edge  of  the  added 
pulse)  is  obtained.  Although  the  maximum  value  of  the  pressure  coefficient  has  decreased 
with  additional  cycles,  the  absolute  maximum  centerline  velocity  increases,  with  the  three- 
cycle  result  showing  a  50-percent  increase  over  that  observed  for  a  single  cycle.  This  result 
implies  that  waveform  can  be  influential  in  determining  the  fluid  response  to  a  pulsation. 

The  radial  velocity  distribution  for  the  three  eases  is  shown  in  Fig.  23  for  a  single  time  of 
the  cycle  (r  =  90  deg).  While  all  the  velocity  profiles  exhibit  a  characteristic  annular  shape, 
significant  variations  in  amplitude  are  obtained  for  differing  repetition  rates.  It  is  noted, 
moreover,  that  as  the  repetition  rate  increases,  the  local  velocity  maximum  occurs  nearer  the 
wall,  indicating  that  the  vibrational  Reynolds  number  is  effectively  higher  than  the  basic 
value  for  the  case. 
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Figure  22.  Centerline  velocity  response  to  asymmetric  waves. 


5.0  CONCLUSIONS  AND  HKCOMMKNDATIONS 

An  analysis  was  performed  and  computer  programs  were  written  to  allow  assessment  of 
the  velocity  field  in  the  fully  developed  region  of  a  pipe  Tor  laminar  or  turbulent,  steady  or 
pulsating  flows.  For  pulsating  flows  no  restriction  on  waveform  or  the  origin  of  the 
disturbance  is  imposed  except  that  it  be  applied  longitudinally. 

For  any  waveform  investigated,  it  was  shown  chat,  in  the  core  of  the  pipe,  the  fluid 
velocity  at  any  time  during  a  cycle  is  an  integral  transform  of  the  imposed  pulsation.  Near 
the  tube  wall,  viscous  effects  become  important,  and  the  velocity  may  exceed  that  observed 
on  the  centerline.  The  extent  of  this  annular  flow  region  is  dictated  solely  by  the  vibrational 
Reynolds  number,  the  effect  being  confined  nearer  the  tube  wall  as  Rev  increases.  The  value 
of  the  maximum  velocity  relative  to  the  core  value  varies  in  a  complex  fashion,  initially 
increasing  with  frequency  but  decreasing  for  high  frequencies.  The  overall  level  of  the 
transient  velocity,  however,  varies  directly  with  a  generalized  amplitude  and  decreases  with 
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Figure  23.  Distribution  of  transient  velocity  with 
radial  location. 

vibrational  Reynolds  number.  Limited  comparisons  of  the  analysis  with  published  data 
obtained  for  a  vibrating  tube  show  good  agreement  for  no  mean  flow  and  laminar  mean 
flow  cases.  In  turbulent  mean  flows,  however,  the  analysis  agreed  with  the  data  at  some 
times  during  a  cycle  but  overpredicted  the  effect  of  pulsations  at  others.  A  crude  estimate  of 
the  effect  of  alteration  of  the  turbulent  exchange  mechanisms  was  found  unsatisfactory  For 
explaining  the  discrepancy. 

A  limited  study  of  the  effect  of  waveform  on  fluid  response  showed  that  the  above 
general  observations  were  valid  for  any  waveform.  For  symmetrical  waves  it  was  found  that 
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the  magnitude  and  radial  or  temporal  variations  in  velocity  were,  in  a  practical  sense, 
virtually  independent  of  waveform.  Asymmetric  waves,  however,  were  shown  to  exhibit 
significantly  different  behavior. 

The  computer  programs  developed  for  the  study  were  found  to  be  satisfactory  although 
several  improvements  arc  warranted.  First,  a  test  of  the  convergence  of  the  Fourier  series 
and  velocity  series  would  be  of  benefit  in  reducing  the  arbitrary  waveform  execution  lime. 
Secondly,  several  features  of  the  harmonic  oscillation  code,  such  as  strength  of  pulsation, 
rms  velocity,  and  a  vibrational  Reynolds  number  calculation,  should  be  included  in  the 
arbitrary  waveform  code.  Finally,  both  programs  should  be  extended  to  provide  computer 
plotting  of  the  results.  It  should  be  noted,  moreover,  that  the  steady  flow  routine 
incorporated  in  either  program  affords  a  simple,  fast,  and  accurate  method  for  computing 
the  velocity  distribution  in  fully  developed  laminar  and  turbulent  flows. 

The  present  study  has  clearly  shown  that  the  complexity  of  time-dependent  flows  for 
even  a  relatively  simple  case  precludes  the  extrapolation  of  steady  flow  results  into  the 
transient  domain.  It  is  recommended,  therefore,  that  work  be  continued  at  AEDC  on 
transient  flow  phenomena.  Two  (ideally,  parallel)  approaches  are  suggested:  (1)  the  present 
“exact”  analysis  should  be  extended  to  (a)  treat  the  developing  pulsating  flow  in  a  tube  and 
(b)  account  for  turbulent  exchange  mechanisms  by  the  method  of  Appendix  A;  (2)  time- 
dependent  numerical  solutions  of  the  governing  equations  should  be  developed  since  they 
offer  the  only  practical  long  range,  general  solutions.  In  addition,  consideration  should  be 
given  to  longer  range  projects  which  (1)  treat  compressibility  phenomena  in  transient  flows 
and  (2)  obtain  basic  data  on  the  effect  of  flow  oscillations  of  various  waveforms  on  the 
velocity  field,  turbulent  structure,  and  skin  friction.  The  first  study  is  of  importance  since 
compressibility  effects  can  have  a  significant  effect  on  the  velocity  field  development 
because  wave  propagation  characteristics  can  be  altered  and  secondary  flows  induced.  An 
experimental  program  is  required,  moreover,  to  serve  as  a  guide  for  modeling  the  effect  of 
pulsations  on  turbulent  structure.  In  particular,  little  or  no  data  exist  on  the  effect  of 
waveform  on  the  turbulent  structure  in  pulsating  flows. 

Finally,  it  has  been  shown  that  How  pulsations  can  have  a  profound  effect  on  heat 
transfer.  Although  myriad  experiments  have  been  performed  and  analyses  presented,  no 
definitive  conclusions  can  yet  be  offered  on  the  effect  of  pulsations  on  heat  transfer  (Refs.  4 
and  8).  An  understanding  of  the  mechanisms  involved  could  possibly  lead  to  revisions  in  the 
design  of  heat  exchange  devices. 
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APPENDIX  A 

FORMULATION  FOR  ARBITRARY  EDDY  VISCOSITY  MODEL 

The  radial  distribution  function  for  the  transient  velocity  component,  f(ij),  is  given  by  an 
equation  of  the  form 


ctf"H  /3f'—  iyf  =  -fi  (A-l) 

where  primes  denote  differentiation  with  respect  to  t}  and 

«(*)  =  1  h  — —  (A-2) 


Pit})  =  i  b(v)  + 


(A-3) 


y  =  Rcv 


(A-4) 


In  general,  f(i))  is  a  complex  function  and  may  be  written 

f(ij>  =  fr(^)  +-  i  f  |(  77) 


(A-5) 


where  fr(?j)  and  f;(7j)  are  the  real  and  imaginary  parts,  respectively,  of  the  function. 

Substituting  Eq.  (A-5)  into  Eq.  (A-l)  and  requiring  a,  a,  (3,  and  7  to  be  real  functions 
leads  to  the  following  simultaneous  equations: 


-b  =  afr"+-  +  yC 


(A-6) 


r 


r 


(A- 7) 


The  real  part  of  the  distribution  is  thus  determined  once  f,  and  its  derivatives  are  known. 

Differentiating  Eq.  (A-7)  and  substituting  in  Eq.  (A-6)  leads  to  the  following  fourth- 
order  differential  equation: 

A<i,)fi(IV)  +  BtJ7)fs<Iin  t  C(ij)  fi(,l)  +  Dl^f;411  *Efj  =  -F  (A-8) 
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where 


A  It))  =  a~ 

BCtj)  =  2a(a'+  f3) 

C (17)  =  a(a~+  2/3'J  h/9(o 
D(t,)  =  aft"  /8j8 ' 


(A-9) 


Thus,  if  €fn|  can  be  described  by  any  analytic  function  with  continuous  second  derivatives, 
the  coefficients  are  uniquely  determined.  Equation  (A-8),  moreover,  may  be  solved  by  any 
of  several  numerical  methods  for  ordinary  differential  equations  (Ref.  22).  The  solution, 
however,  will  require  initial  estimates  of  f,  and  its  derivatives  at  one  or  more  points  in  the 
flow  and  satisfaction  of  the  boundary  condition  fj(l)  =  0.  These  estimates  may  readily  be 
obtained  from  the  em|  =  constant  solution  of  Eq.  (80). 

One  solution  procedure  which  has  been  programmed  in  Fortran  IV  involves  reducing  the 
order  of  Eq.  (A-8)  by  introducing  the  relations 

with 


Equation  (A-8)  is  thus  reduced  to  a  series  of  first-order  simultaneous  equations: 


■  4  [cs3,  c:e2 *  d8j  .  Et. .  f] 


(A- 11) 
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*3 


3  p2 


l' 


(A-ll) 

Concl. 


Integration  of  each  equation  was  accomplished  using  a  fourth-order  Runge-Kutta  method  in 
conjunction  with  the  secant  method  for  simultaneous  nonlinear  equations  (Ref.  23).  The 
real  part  of  the  distribution  function  is  then 


f 


r 


a  .  £ 

y  g2 T  y  gi 


(A- 12) 


The  transient  velocity  at  each  point  of  the  flow  field  is 

Vj(V)  =  1^(37)  cns  r—  Mtj)  sin  r  (A-13) 

and  the  instantaneous  velocity  profile  is 

V r)  =  -  Vj  (77, r)  (A- 14) 


with  V0(ij)  given  by  Eq.  (46). 
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APPENDIX  B 

A  NOTE  ON  COMPLEX  FOLRIER  EXPANSIONS 

The  Fourier  trigonometric  expansion  of  an  arbitrary  periodic  function,  f(t),  which 
satisfies  the  Dirichlet  conditions  (Ref.  24)  may  be  written 

A|>  SO  « 

f  (t>  =  —  ^  2]  A  cos  neul  t-  V  lin  sin  n&jl  (B-l ) 

2  “  1  "  n“ 


where  the  Fourier  coefficients  are  given  by  the  transforms 


and 


A 

n 


f(t ')  cos  {ncdl ')  dt';  n  -  0.  1.  2.  . . . 


(B-2) 


13 

II 


f  (l  ')  sin  (nojt')  dt  n 


]  2 

J  f  **f  •  *  * 


(B-3) 


The  basic  period  of  the  function  is  P  =  2x/w,  and  ijcj  represent  the  various  harmonics  of 

to. 


it  is  often  convenient  to  express  Hq.  (B-l)  in  its  equivalent  complex  form. 


f(0 


mail 


n  =  *-ac 


(B-4) 


which  is  readily  obtained  by  a  direct  substitution  of  the  complex  definitions  of  the  sine  and 
cosine  into  the  trigonometric  series.  When  that  is  done,  the  coefficients  are  defined  by 


or 


f(t') 


—  i  naj  l  * 
d 


dt' 


(B-5) 


M 


-  5B 


n 


} 


(B-6) 
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The  functions  f(t)  considered  in  ihis  report  are  always  real;  thus  it  is  possible  to  rewrite 
Eq.  (B-4)  in  a  form  that  is  more  amenable  to  computation.  For  this  purpose,  the  complex 
expansion  is  rewritten  as 


fft)  =  jt  -Ik  ciI,fcl  + 

0  -i  n 


oc 


where  Hn  is  still  given  by  Eq.  (B-5)  or  (B-6)  but 


i.iil&n 

e 


(B-7) 


(B-8) 


and 


(B-9) 


Comparing  Eqs.  (B-6)  and  (B-9),  one  sees  that  H„  and  H  _  |  n|  are  complex  conjugates.  The 
terms  to  be  summed  in  Eq.  (B-7)  are,  likewise,  complex  conjugates.  Thus,  for  any  real 
function  f(l),  Eq.  (B-7)  becomes 


f(i)  -  K 


I  K 

n  =  ] 


■  n  d>  1 


(B-10) 


where  it  is  understood  that  n  is  a  positive  integer  and  the  coefficients  are  defined  by  Eq. 
(B-6). 
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A 

An 

Ap 

a* 

Cl 

Cl* 

b 

C 

Cf 

FnO?,t) 

g-(t) 

t'(v) 

f,(v) 

CAv) 

G* 

I,(i|/2x) 

K„(i''2x) 

M,.(x) 

m 

tti 


NOMENCLATURE 

Dimensionless  vibralion  amplitude,  Z0/R 
Fourier  Cosine  Transform,  Eq.  (71) 
Dimensionless  pressure  amplitude,  J]|/Qo 
Frequency  parameter,  (Rev)t/2 
Generalized  amplitude,  Eq.  (40) 

Modified  amplitude,  Eq.  (54) 

Fourier  Sine  Transform,  Eq.  (72) 

Mean  value  of  arbitrary  parameter 
Constant 

Fanning  friction  factor,  u2 
Arbitrary  function,  Eq.  (28) 

Arbitrary  forcing  function,  Eq.  (40) 

Radial  distribution  function,  Eq.  (57) 

Imaginary  part  of  distribution  function,  Eq.  (80) 
Real  part  of  distribution  function,  Eq.  (79) 
Dimensionless  mass  How  rate,  Eq.  (33) 

Complex  Fourier  Coefficient,  Eq.  (62) 

Modified  Bessel  function  of  first  kind 
Modified  Bessel  function  of  second  kind 
Modulus  of  polar  representation  of  l„(il,'2x) 
Empirical  coefficient,  Eq.  (49) 

Mass  flow  rate  (M/t) 
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P 

P 

P' 

P* 

R 

Re 

Re* 

Ren 

Rev 

Re? 

RujUj 

M,z) 
(r',0' ,z') 


nth  harmonic  or  term 

Instantaneous  pressure  (m/Lt2)  or  period  (t) 

Mean  pressure  (M/Lt2) 

Fluctuating  pressure  (M/Lt2) 

Reference  pressure  (M/Lt2) 

Tube  radius  (L) 

Mean  throughflow  Reynolds  number,  2<U>R/n 
Friction  Reynolds  number,  U*R/v 
Harmonic  Reynolds  number,  wZ0R/r> 

Vibrational  Reynolds  number,  uR2/v 

Modified  vibrational  Reynolds  number,  Rev/[1  +  (eni,/j')] 

Correlation  coefficient,  u,'uj7U*2 

Noninertial  cylindrical  coordinate  system 

Inertial  cylindrical  coordinate  system 


S  Pulsation  strength,  Eq.  (82) 

s  Empirical  coefficient,  Eq.  (48) 

T  Time  interval  (t) 

t  Time  (t) 

(U,V,W)  Reference  frame  velocity  (L/t) 

Uh  Harmonic  velocity,  uZ0  (L/t) 

U*  Friction  velocity,  \f < <rw”> /p 

(u,v,w)  Mean  velocity  components  (L/t) 

(u '  ,v '  ,w ' )  Fluctuating  velocity  components  (L/t) 
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V, 

Vc 

Vi 

Vi 

Vo 


V; 


Vi' 

Xi 


Xi.Xj 


X 

<x> 

X 

Zo 

A 

8 

em0 

[ 

r 

V 

An 

M 

Mo 


Dimensionless  transient  velocity 
Centerline  velocity,  W(o,t)/U* 

Instantaneous  velocity,  w(r,t)/U* 

Reference  frame  velocity  component  (L/l) 
Dimensionless  mean  velocity 
Mean  velocity  component  (L/t) 

Fluctuating  velocity  component  (L/t) 

Reynolds  stress  component 
Body  force  (ML/t2) 

Cartesian  coordinates 

T 

Time  mean  value,  1/T  Jq  xdt 
Mean  value  over  a  cycle,  1/2tt  J^xdl 
Spatial  mean  value,  X17  dij 
Vibration  amplitude  (L) 

Cubical  dilitation, 

Boundary-layer  thickness  (L) 

Eddy  viscosity  of  mean  flow  (LVt) 

Eddy  viscosity  of  transient  flow  (L2/t) 
Dimensionless  axial  coordinate,  Z/R 
Dimensionless  radial  coordinate,  r/R 
Fourier  coefficient  relations,  Eqs.  (74)  and  (75) 
Absolute  viscosity,  (M/Lt) 

Bulk  viscosity,  (M/Lt) 
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Kinematic  viscosity  (LM)  or  order  of  Bessel  function 
Phase  of  polar  representation  of  I„  (i1/2x) 

Density  (M/L3) 

Wall  shear  stress,  (M/Lt2) 

Dimensionless  time,  a>t 
Pressure  coefficient,  P  -  P*/pU*2 
Mean  flow  pressure  coefficient 
Transient  component  of  pressure  coefficient 
Frequency  (1/t) 


